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Abstract 

It is shown that the theoretical arguments for the unavoidable collapse of 
massive stars, which are valid within standard general relativity, are not valid 
in extended gravity. In particular a calculation is reported of neutron stars 
within a theory obtained from an Einstein-Hilbert action involving a function 
of — 1/2R^,,R^'^ added to the Ricci scalar R. The calculation suggests 
that within that theory there are stable neutron stars with arbitrarily high 
baryon number. 
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I. Introduction 

Claims for the existence of black holes rest upon observational evidence 
combined with theoretical arguments. Here I analyze whether the theoreti- 
cal arguments, valid within general relativity (GR), remain valid when GR 
is modified, for instance by adding other terms to the Ricci scalar in the 
Einstein-Hilbert action. 

Black holes have a long and interesting history.l^ The realization of a 
Schwarzschild singularity was taken as a serious possibility in 1930, after the 
work of Chandrasekhar. He proved that white dwarf stars of high enough 
mass (and low enough temperature) are unstable against gravitational col- 
lapse. As is known today those stars do not collapse directly to black holes 
but suffer a supernova explosion, a core remaining in the form of a neutron 
star which is believed to evolve towards a black hole if the mass is large. 
The idea of collapsing stars had a strong oposition on the part of leading 
astrophysicists of the time, in particular Eddington who assumed that some 
effect, not yet known, would prevent the collapse. 

The next important date was 1939 when Einstein, then in Princeton, 
and Oppenheimer in California wrote three important papers which I shall 
comment in the following. All three papers dealt with spherically symmetric 
bodies, which might be studied using standard coordinates with metric 

ds^ = exp {a (r)) dr^ + r^dQ^ - exp (/3 (r)) dt\ (1) 

In these coordinates the GR hydrostatic equilibrium equation is 

dp q — P k (m + Attt^p) , ^ /"'" 9 , 

-f + 2^— ^ + -^^— — ^ p + p =0,m= / Anr^pdr, 2 

dr r oTrr^ — Zkmr Jo 

where k is Sir times the Newton constant, r the radial coordinate, and p,p 
and q are the mass density, the radial pressure and the transverse pressure, 
respectively. Oppenheimer and Volkoff (OV)^ solved eq.(|2]) with the assump- 
tion of local isotropy, that is p = q, and the equation of state, p = p{p), 
appropriate for a fluid of free (non-interacting) neutrons. The result was 
an one-parameter family of equilibrium configurations where the total mass, 
M = m {00) , increases with the central density until a maximum about 0.7 
solar masses (the OV limit) and then decreases. The relevant conclusion was 
that there are no equilibrium configurations of free neutrons with mass above 
the OV limit. Later on it has been proved that for interacting neutrons a 
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limit still exists although it is somewhat higher, about 2 solar masses. Thus 
the question was, what happens to the core of a white dwarf after a super- 
nova explosion when the core has a mass above the OV limit?. In order to 
answer this question Oppenheimer and Snydei^ studied the evolution of a 
spherical cloud of dust. The result was that whatever is the initial config- 
uration the whole cloud collapses to a black hole. This suggests that any 
neutron star with a mass above the OV limit would suffer a collapse, quali- 
tatively similar to the cloud of dust. Later calculations have confirmed this 
assumption, which has lead to the current wisdom that all observed compact 
(cold) objects with mass above the OV limit are black holes.l^ 

EinsteirP studied a model of star consisting of a spherical distribution of 
particles each one moving in a circle around a central point. Thus the matter 
in the star may be considered as a fluid with extreme local anisotropy, that is 
nil radial pressure, p = 0. If this is put in the hydrostatic equilibrium eq.(l2]) 
we get 

, , km (r) . , . , 

^W = ^ T^r^pir). 3 

ovrr — 2km [r) 

which has a solution for any (spherical) density distribution. In fact, given 
any function p (r) we may get m(r) via the second eq.(I2]) whence eq.(|n]) gives 
q (r) . Actually those functions p (r) which lead to a pressure violating the 
constraint g (r) < p (r) should be excluded because this would mean that 
some of the constituent particles move with a velocity greater than that of 
light, c. On the other hand it is easy to prove that the solutions of eq.(l3]) 
with g (r) < p (r) are stable against homologous collapse, which may be 
represented by a family of configurations depending on a parameter A such 
that 

r = Ax, p (Ax) = X^'^p (x) , q (Ax) = \~^p (x) , m (Ax) = m (x) . 
That is we should consider a sequence of configurations fulfilling 

q (Ax) q (x) km (x) 

p (Ax) p{x) Stt Ax — 2A;m (x) ' 

with decreasing values of A. It is obvious that the ratio q/p will increase 
indefinitely for A — )■ km{x) /{Attx). Then the homologous collapse should 
stop before the constraint q{r) < p (r) is violated. From this result Einstein 
concluded that the existence of a limiting velocity, c, prevents the collapse 
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to a Schwarzschild singularity, and he even asserted that for any star model 
(i. e. not necessarily fulfilling p = 0) this would be the caseP 

As is well known today Einstein was wrong, because in the real world 
particles interact with each other and the interactions tend to produce local 
isotropy. Actually between local isotropy, p = q, assumed by OV^ in their 
neutron star calculation and extremal local anisotropy, p = 0, as in Einstein 
model, there are intermediate cases. Indeed calculations have been performed 
which prove that for some amount of local anisotropy, q ^ p, the OV limit 
may be surpassed in neutron starsP However it is easy to prove that stars in 
stable equilibrium should have local isotropy, which in particular shows that 
Einstein's star model is unstable. In fact, let us assume an equilibrium con- 
figuration with density distribution p (r) , p{r) and q{r) being the radial and 
transverse pressures respectively. Now, in the absence of external electric or 
magnetic fields a fiuid of neutrons with given baryon number density, n (r) , 
has the minimal energy density for an isotropic velocity distribution, that is 
p{r) = q{r). As a consequence if we start with an equilibrium star config- 
uration presenting local anisotropy and change it to another configuration 
with the same n (r) but local isotropy everywhere, then the energy density 
decreases at every point where previously q ^ p. In summary we pass from 
an equilibrium configuration with local anisotropy to another configuration, 
maybe out of equilibrium, with the same baryon number but a smaller total 
mass. Consequently the equilibrium star configuration is not stable. 

Many authors have been reluctant to accept the existence of a spacetime 
singularity in the line of Eddington and Einstein. There are good reasons 
for that. In fact according to GR the metric of spacetime is ruled by second 
order differential equations and therefore it is natural to ask that the metric 
coefficients are twice derivable everywhere, with no singularity. The current 
solution to the problem is to assume that the singularity which appears in 
classical GR may be cured by quantum effects. But the good solution would 
be that quantum effects could modify GR in such a way that stars never 
collapse. On this line there is a recent proposal of "gravastars" , that is stars 
with a core where the vacuum, with equation of state p = —p, provides a 
strong repulsive pressure able to stop the collapse.'^ In this paper I propose 
a more natural solution, namely that the correct gravitational theory is an 
extended GR, derived from a modified Einstein- Hilbert action,'^ and that in 
such gravity theory stars never collapse. The purpose of this paper is to 
provide arguments suggesting that this possibility is worth to be seriously 
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studied. 



II. Slightly relativistic stars in f(R) gravity 

The most popular stars in Newtonian gravity are polytropes, that is stars 
consisting of a fluid with equation of state 

p = Kp^, 

where K and 7 are constant. It is well known that polytropes are stable 
(unstable) if 7 is greater (smaller) than 4/3. It is the case that the most im- 
portant types of Newtonian (or slightly relativistic) stars where gravitational 
collapse is currently assumed are (massive cold) white dwarfs and superma- 
sive StarBand both may be treated as close to polytropes with 7 = 4/3. Thus 
these stars are at the limit of stability and it is believed that they become 
unstable due to GR corrections. In fact for a slightly relativistic star eq.(I2]) 
may be written, assuming local isotropy, that is p = q, 

dp k (m + Airr^p) , , kmp , Anr^p p km , , , 

T = - 02 ot. iP + P) - -^(1 + + - + (4 

ar ovrr^ — Zkmr STrr^ m p Airr 

The main term in the right side correspond to the Newtonian approximation 
and the remaining 3 terms are the GR corrections. It is easy to see that every 
GR correction contributes to increase the gravitational field with respect to 
the Newtonian one, therefore inducing unstability towards collapse. In the 
following I show that corrections derived from f{R) gravity have precisely 
the opposite effect. 

F{R) gravity derives from the Einstein-Hilbert action 

S=jj^J d^x^g {R + F) + Smat, (5) 

where R is the Ricci scalar and F (R) is an arbitrary function or R, which I 
shall assume twice derivable. The action eq.(l5]) gives the field equatioiJ^ 

(1 + F^ {R))R^, - ^g^u (R + F (R)) = V^V.Fi - g^.F^ + kT^,, (6) 

where Fi = dF/dR and R^u (T^u) is the Ricci (stress-energy) tensor. The 
divergence of eq.(l6]) leads to 

kWT^, = R,.WF^ (R) + (V, - V,) Fi {R) . (7) 
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For a static spherical body in the metric eq.(IT]) only the component with 
u = 1 of the tensor eq.(I7]) gives a nontrivial result, that is (for local isotropy 
p = q) 

% + \P\P + P) = k-'F, (R) R' + 1/3" - - ^) , (8) 

where F2 = d?F/dR^,R' = dR/dr,P' = dp/dr,(3" = d^(3/dr^,a" = d^a/dr\ 
Now let us assume that the deviation from GR is not large, that is the right 
side of eq.® is in some sense small (note that the left side equated to zero 
just gives the GR eq.(l2])). Then for a slightly relativistic star it is consistent 
to neglect the terms a" and (3" in the right hand side and we get 

I + l^'iP + P)^ -\P\P + P) + k"'F, (R) R' [rI - ^) . 

Similarly we have, from Einstein's equation, 

Rl = iA;(p-p) ~ ^kp, 

where I have taken into account that p << p for an almost Newtonian star. 
But in addition we have 

4 o km 2 
-npr < m, — « — , 

the former inequality being a consequence of the fact that the density de- 
creases with increasing radius and the latter is true for slightly relativistic 
stars where km « Airr. Thus we may neglect the term with Rl whence, 
taking eq.Q into account eq.([H]) , gives finally the following effective gravi- 
tational field 

p dr ovrr^ m p Airr r'^p dr 

The relevant result is that the latter term, which is the correction derived 
from f{R) gravity, is always positive because dp/dr < and the inequality 
F2 (p) > holds true for any viable f{R) theory .f^^ In conclusion any modifi- 
cation of GR via f(R) theory tends to compensate for the unstability created 
by GR in slightly relativistic stars. Whether this is enough to prevent the 
gravitational collapse of white dwarfs or supermassive stars may depend on 
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the form and size of the function F{R), but I will not pursue the study here. 
In summary, although I cannot claim that f{R) theory always prevents the 
gravitational collapse of white dwarfs and supermassive stars, it is true that 
their collapse is not unavoidable in f{R) theory, at a difference with standard 
GR. 

III. Neutron stars in quartic gravity 

Our second example will be the case of neutron stars in quartic gravity, 
where the Lagrangian F of eq.® i^ 

F = aR^ + bRf.^R^"', a = -6/2 = Ikm^. (10) 

The gravity theory resulting from this choice is not viable for, at least, two 
reasons. Firstly it contradicts Solar System and terrestrial observations or 
experiments, which require y/a and yZ—b not greater than a few millimeters. 
Secondly the weak field limit of the theory presents ghosts. The difficulties 
may be solved if we assume that the theory resting upon eqs.(l5]) and ffTOl) is 
an approximation, valid for fields of the order of those appearing inside or 
near neutron stars, of a theory giving a negligible correction to GR for weak 
fields. For instance, we might assume that eq.l JTOj) is an approximation of 
the theory resting upon the Lagrangian 

F = aR^ + hR^^R^"" - c log (l + {a/c)R^ + {h/ c)R^,R^"') . (11) 
In particular if we choose c ~ 1/(10^° km^) eq.ffTTl) becomes 

F ~ ^ (aR^ + bR^uR^^y ~ 10-^^ R for the Solar System, 

F ~ aR^ + bRf.^R'"' ~ lO'^R for neutron stars, (12) 

and the relative error due to the terms neglected in going from eq. (fTT|) to 
eq.( IT2l) is smaller than 10~^^ in both cases. I have taken into account that 
R^ ~ R^yR^^ ~ {kpY and that typical densities are p ~ lO'* kg/m^ in the 
Solar System and p ~ 10^^ kg/m^ in neutron stars. I conclude that, in this 
example, the correction to GR would be negligible for the Solar System (and 
even more so for the cosmology at present time), but quite important for 
neutron stars. In any case it is irrelevant for the purposes of this paper 
the question whether the theory is viable or not. Indeed I do not pretend 
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to study neutron stars with the "correct" gravity theory which is unknown, 
but to provide a counterexample to the conjecture that stable neutron stars 
cannot exist with more than twice the baryon number of the Sun. The coun- 
terexample proves that such a limit does not hold true in extended gravity 
theories. 

In the following I report on the results obtained from the calculation of 
noninteracting neutron stars using the gravity theory defined in eq. ( ITOj) . The 
field equation i^ 



G^u + a 



-GG^iy + -^g^u{G^ — GxaG'^") — G'^Ga 



+2" 



G 



V.V.G^ - V.V^G: + G^ 



kT, 



(13) 



written in terms of the Einstein tensor G^,^ = Rfj^^ — \g^uR rather then 
the Ricci tensor Rav When the components of Guy are written in terms 



of the parameters a and /3 of the metric and their derivatives^ we get a 
system of 3 independent fourth order differential equations (one of the four 
components of the tensor eq.f fT3|) is not independent because in spherical 
symmetry G22 = G33.) The constraints that a and /3 are finite at the origin 
and go to zero at infinity fix all the initial or boundary conditions except 
one, so that the solution of the 3 (coupled, nonlinear) ordinary differential 
equations provides a one-parameter family depending of the central value of 
the component G\. For simplicity in the calculation I have used as equation 
of state the following 

Pmat = ^Pmat + GpTat, G = 8.17 X 10^ kg^/^m-^/S, (M) 

which has the correct behaviour for a fluid of free neutrons both at high 
and low densities. The constant G is so chosen that a general relativistic 
calculation gives the same result as the original Oppenheimer and VolkofP 
for the maximum mass stable star (see Table 2 below). 

A relevant quantity is the baryon number of the star, A^, which may be 
calculated from the baryon number density n{r) via 



N 



n{r)r'^dr 



^1 - km{r) / (Anr) 



(15) 



The number density n may be related to the mass density and the pressure 
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from the solution of the equation 



dpmat 

Pmat ^ ; Pmat) 

an 

which follows from the first law of thermodynamics. Inserting eq.( !T^ here 
we get a differential equation which may be easily solved with the condition 
pmat In — !■ for p — )■ 0, /i being the neutron mass. I get 

n = C'''vTat{^vTat + cf". (16) 

The details of the calculation may be seen elsewhere.'^ In Table 1 I give 
the results obtained. For comparison in Table 2 the results of a calculation 
within GR are presented. 

Table 1. Calculation within extended gravity theory. The parameter k is 
Stt times Newton's constant. The central matter density, Prn.at(0), is in units 
Pc = 7.2 X 10^^ k.g/w? and similar units are used for the products of k~^ 
times the components of the Einstein tensor, GJ^, at the star center. Central 
baryon number density, ?t.(0), is in units pd p, p being the neutron mass. The 
star radius, R, is in kilometers, the mass, M, in solar masses and the baryon 
number, A^, in units of solar baryon number. I report also the dimensionless 
fractional surface red shift, A A/ A = 1/ ^Jl — 2M/R — 1, and the percent 
binding energy, BE = 100(A^ — M)/N. An expressions like 1.QE2 means 
1.6 X 10^ 



k-^G\ (0) 
k-'Gi (0) 


0.01 


0.1 


1 


10 


100 


1000 


10000 


0.18 


0.82 


4.5 


34 


3.1E2 


3.0E3 


3.0E5 


Pmat (0) 


1.6E2 


2.5E3 


4.3E4 


7.8E5 


1.QE7 


3.2E8 


6.3E9 


n(0) 


56 


4.5^2 


3.7 E3 


3.3^4 


3.2E5 


3.0EQ 


2.8E7 


R 


10.7 


6.7 


4.0 


2.7 


2.1 


2.2 


2.2 


M 


0.67 


0.80 


0.60 


0.39 


0.264 


0.268 


0.292 


N 


0.73 


0.94 


1.03 


1.13 


1.44 


2.00 


2.63 


BE 


8.9% 


15% 


41% 


65% 


82% 


87% 


89% 


AA/A 


0.106 


0.22 


0.34 


0.31 


0.23 


0.23 


0.26 



Table 2. Calculation within general relativity. The meaning of all sym- 
bols is as in Table 1, but here the matter pressure and density equal k~^ 
times the appropriate components of Einstein's tensor. 
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k-^Gl (0) = p{0) 


0.01 


0.04 


0.1 


1 


10 


100 


1000 


0.18 


0.46 


0.89 


5.3 


39 


3.4 X 10^ 


3.1 X 10^ 


R 


11.9 


9.6 


8.3 


5.8 


5.2 


6.6 


6.6 


M 


0.67 


0.72 


0.70 


0.55 


0.39 


0.40 


0.44 


N 


0.69 


0.74 


0.73 


0.55 


0.36 


0.37 


0.42 


BE 


2.9% 


3.4% 


3.4% 


-0.73% 


-8.1% 


-8.0% 


-5.7% 


AA/A 


0.094 


0.13 


0.15 


0.18 


0.13 


0.11 


0.12 



IV. Discusion 

The most dramatic difference between the calculations within extended 
gravity and within GR is in the variation of the baryon number with the cen- 
tral density. In GR there is a maximum about = 0.74 times the baryon 
number of the Sun for some central density, p (0) , and after that it decreases 
with increasing p(0) (which gives the OV limit for the baryon number). In 
our calculation with modified gravity N increases monotonically with the 
central density, which suggets that there is no limit. Of course for the ex- 
tremely high densities which appear in Table 1 the individual neutrons would 
not exist and a fluid of different particles would appear, but we have used 
the same equation of state (of free neutrons) throughout. It is unlikely that 
a change of the equation of state could produce a very different qualitative 
behaviour, taking into account that for high density the equation should have 
the form p ^ p/3 in any case. The binding energy is also dramatically differ- 
ent in both calculations. Finally the gravitational surface redshift surpases 
0.3 whilst in GR it is always smaller than 0.2, which seems to agree better 
with observations.'^ The comparison with observations is however not too 
relevant because the results do not take into account the interaction between 
neutrons. 

If the "correct" gravity theory gives results qualitatively similar to those 
in Table 1, the evolution of any massive star would end in a compact star 
with relatively small mass but a very large baryon number. The essential 
conclusion of the paper is that the theoretical arguments supporting the 
existence of black holes, which rest upon standard GR, may be no longer 
valid in extended gravity theories. If it turns out that GR must be modified 
(possibly by quantum effects), then it might be the case that the alledged 
observational evidence for black holes should be revised. 
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